A proposal for measuring correlation functions in interacting gases 
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We propose an experimental method that allows for scaling of all the diagonal correlation functions 
(such as the density and pair correlation function) during the expansion of harmonically trapped 
cold atomic gases with arbitrary interactions even if the trapping is anisotropic. To do this we 
require that the interaction be tunable by the use of a Feshbach resonance and that the trapping 
potential can be manipulated independently along the axial and radial directions. This permits a 
much greater accuracy and resolution in measurements of correlation functions. We explain the 
explicit control procedure, outline the scaling ansatz that we use and briefly discuss the systems 
where it could be applied as well as possible experimental obstacles to its implementation. 

PACS numbers: 



How can we measure correlation functions in harmon- 
ically trapped cold atomic gases? In order to measure 
them with good accuracy, the trapped gas is first al- 
lowed to expand ballistically to improve optical resolu- 
tion and reduce its optical depth. However the inter- 
actions present in the gas can in principle modify the 
correlation functions significantly during the expansion 
and it is far from obvious how to relate the measured 
values after expansion to those that existed in the trap 
before the expansion [l| . Here we propose a new exper- 
imental method whereby the release from the trap of an 
interacting gas is accomplished in such a way that cor- 
relation functions exhibit scaling during the expansion, 
and therefore their initial values can be measured with 
great accuracy. We shall not discuss the off-diagonal cor- 
relation functions such as the one-body density matrix 
(■0'l'(r)V'(r')) with r ^ r'. However the present method is 
applicable to all real space diagonal correlation functions. 

In most experiments that measure the density n(r), 
the atomic gas is allowed to expand and, when it is large 
enough for a satisfactory optical resolution, it is imaged 
destructively. In many cases (e.g. for ideal gases or in- 
teracting ones in the Thomas- Fermi limit), the density 
after expansion is then related to the one before expan- 
sion using scaling transformations [3, Q . By scaling we 
mean a coordinate transformation of the kind: 



\x{t)x, y \{t)y, z Xz{t)z 



(1) 



such that the density at a time t after release from the 
trap can be written in terms of the density at t = 0: 

n{r,t) = ano{x/Xx,y/Xy,z/Xz) (2) 

The Xi(t) are known as scaling functions and a = 
{XxXyXz)~^ is a time dependent normalization factor. 



While the density gives us very important information 
about such things as thermodynamic properties and col- 
lective oscillations, it does not reveal all the many-body 
aspects of the system. In fact we can probe new physics 
by looking at other quantities such as the pair correlation 
function which is defined as: 



n(ri,<)n(r2,i) 



(3) 



with 
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n(^)(ri,r2,t) = {^^ {v,,t)^^ {v2,t)^{v2,t)i,{v,,t)) (4) 

where the ?/; are atomic field operators. The pair corre- 
lation function g(ri,r2) is related to the probability of 
finding an atom at T2 given that there is an atom at ri. 
If there is more than one spin species present (as in the 
case of the superfiuid Fermi gas) then we can consider 
also the opposite spin pair correlation function (see e.g. 
jig] ) where the atoms at ri and Y2 are in different spin 
states. The pair correlation function has been the subject 
of recent experimental studies [3, [1, [1, S while others 
have focussed on correlations in momentum space [8, ^ . 
The pair correlation function can reveal new physics that 
is not measured by the density. For example, it can re- 
veal the "bunching" of thermal bosons Ji] , a phenomenon 
which is also predicted to occur in the unitary Fermi su- 
perfiuid [l0|. The difficulty with measuring g^'^\vi,V2) is 
that it has important features at small values of |ri — r2|, 
sometimes comparable to the interatomic distance. This 
requires very good optical resolution which usually can 
only be achieved through an expansion of the gas. 

This leads us naturally to the following question: as far 
as real-space correlation functions arc concerned (such as 
the density or the pair correlation function) , to what ex- 
tent can they be faithfully preserved during the expan- 
sion in such a way that a subsequent measurement will be 
able to recover the values that existed in the trap prior to 
the expansion? For example, if we measure g^^-' (ri, r2, t) 
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after an expansion time t, can we reconstruct the initial 
g^^-* (ri, r2, 0) before the release from the trap? 

Ideally we would like to have a scaling relation like 
that of Eq.([21) for g'^-^\ Then we would have simply that 
(jf^^^(ri,r2,i) oc g^^^rii/Xi,r2i/\i,0). Could such a scal- 
ing relation exist? Unfortunately the answer is no in 
general, as shown in p^ : in an atomic gas, the exis- 
tence of interactions imposes certain boundary conditions 
on the many-body wave function and, as we shall now 
show, these conditions are incompatible with scaling for 
all correlation functions except the density. At distances 
shorter than the interatomic separation, the wave func- 
tion is determined up to a proportionality constant by 
the Bethe-Peierls boundary condition through which the 
interactions appear: for any value of the scattering length 
a the many-body wave function \[' obeys the condition for 
all pairs of atoms ij: 



a r,, - 



(5) 



where the limit is taken keeping all other atoms and the 
center of mass -|- rj fixed. Having taken care of the in- 
teractions via the boundary conditions, the only other 
condition on "if is that it must obey the free-particle 
Schrodinger equation [Tl| . These boundary conditions 
apply also to all diagonal correlation functions: for ex- 
ample (7(^)(ri,r2) obeys the same condition as a func- 
tion of |ri — r2|. The only exception is the density n(r) 
since we need at least two coordinates to probe what hap- 
pens as two atoms come close to each other whereas the 
density depends on a single coordinate. Eq.([5]) has two 
characteristic features: I) it is isotropic (since it depends 
on the modulus of — rj and not on its direction) and 
II) it introduces a length scale a. These two features will 
pose problems for any scaling solution of the form Eq. ([T]) : 
I) requires that all scaling functions Xi{t) be identical, 
which is in general incompatible with the free expansion 
from an anisotropic trap. II) means that, even if the trap 
is isotropic (leading to identical Xi{t)), a rescaling of the 
coordinates would change the ratio a/\ri — rj\ which is 
required to be fixed. These problems do not appear in 
the ideal gas case where a = 0. If a — > oo (e.g. the case 
of the unitary Fermi gas), II) is of no consequence but I) 
still requires that the trap be spherical. 

Nevertheless we can show that there exists in prin- 
ciple an experimental procedure that can restore scal- 
ing during expansion to the many-body wave func- 
tion of a gas initially in equilibrium, trapped in an 
anisotropic trap and interacting via an arbitrary scatter- 
ing length a. By this we mean that |^(ri, rjv, t)P oc 
|*(ri/A(i), ...,rA,/A(i),0)p. Note that information 
about the initial off-diagonal properties is lost and so 
we cannot for example access the momentum space cor- 
relation functions in the same way using this method. 

To defeat problem I) we notice that it is not the initial 
trap anisotropy (i.e. the difference between the Ui) which 
destroys scaling but rather the expansion anisotropy (the 
difference between the Xi{t)). A judicious time depen- 



dent modulation of the trapping constants iOi{t) during 
the release process is sufficient to restore isotropy (all di- 
rections scaling with the same X{t)) as we shall see below. 
Problem II) can also be dealt with if the scattering length 
can be changed during the expansion by modifying the 
applied magnetic field - which requires the existence of 
an accessible Feshbach resonance. The idea is to change 
the scattering length a as a function of time during the 
expansion so as to preserve Eq. ([5]) . 

We now briefly discuss the calculations which underlie 
the scaling approach. Indeed, they follow almost exactly 
those in [l^. The main differences with respect to that 
Reference are: here we consider a general anisotropic 
trap (as opposed to a spherical one), a different time- 
dependence of the trapping frequencies uJi (t) and an ar- 
bitrary initial scattering length, instead of the unitary 
limit. Therefore we shall give here only an outline refer- 
ring the reader for any further details to the more com- 
plete discussion in that Reference. 

For concreteness let us specialize to the experimentally 
relevant case of an initially cigar-shaped potential (even 
though this assumption is not necessary for our proce- 
dure since we could envision other possibilities such as a 
pancake geometry): 

U = im^i(t) {x^ + y^) + ^mul{t)z'' (6) 

where for t < we have and LOz < uj±. 

We also assume that, prior to expansion, the gas is in a 
(many-body) eigenstate \I>(ri, tat, i = 0) in the trap 
U{t — 0), obeying the boundary conditions Eq.® for a 
given scattering length a. We make an ansatz for the 
time-dependence of the wave function: 

/ ^ 

^'(ri, ...,rN,t) = 7V(<)cxp i'^mrjX/2hX 

x*(ri/A,...,rAr/A,t = 0) (7) 

where there is a time-dependent normalisation M and a 
single scaling factor X{t) which is the same for all di- 
rections. For a general value of A(<) this wave function 
no longer respects the boundary conditions for the ini- 
tial value of the scattering length and so we must change 
them using a Feshbach resonance according to the ex- 
pression 



a{t) = a(0)A(i) 



(8) 



Since A(t) generally increases with time, the scattering 
length must also increase (in modulus) but the dilu- 
tion parameter n\a\^ (where n is the density of the gas) 
will remain constant. This procedure solves problem II). 
Now we insert this wave function into the free-particle 
Schrodinger equation, use the fact that ^{t = 0) is an 
eigenstate of the harmonic potential and, by equating 
separately terms in z^^'(i = 0) and in x'^'i'{t = 0), we 
find that it is indeed a solution provided that we can 
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satisiy the two equations 



A = 



A 



_(0) 



A3 

A3 



(9) 



If a physically acceptable solution X{t) of these equa- 
tions can be found then also difficulty I) will have been 
overcome. Indeed such a solution does exist if the trap- 
ping constants ujz.i_{t) can he changed independently in 
an appropriate manner. For example, one such solu- 
tion is to simply switch off the trapping potential along 
z instantaneously as usual {LOz{t > 0) = 0) implying 
\{t) = ^\ + tj40)2i2 (with A(0) = 1 and A(0) = 0) 
which in turn requires that 



c^±(0) 

Vc^i(0)-t^^(0) 



t < 



t > 



(10) 



i.e., an initial discontinuous change of w_l followed by a 
gradual decrease to zero going as 1/i^ at large times. It is 
clear that, for highly anisotropic traps, the discontinuity 
can be neglected (see Fig. [1}. 

So, to summarize, the experimental procedure con- 
sists of simultaneously changing during the expansion 
the radial trapping potential according to Eg. pH]) and 
the scattering length according to Eq.®, i.e. a{t) — 
a(0)v^l+a;,(0)2i2. 

We see therefore that the restriction due to the bound- 
ary condition is not on the initial anisotropy of the trap 
but on the anisotropy of the expansion. Here it pro- 
ceeds isotropically for all times and |^'(i)P expands in a 
self-similar way so that the expansion acts as a perfect 
microscope, enlarging without distortion all the features 
of the cloud. We propose therefore that this method, 
when applicable, be considered the "best practice" for 
measuring real-space correlation functions. 

To illustrate a practical implementation, we consider 
the example of a Fermi gas of atoms, in a mixture of 
two different spin states. We consider an optical trap, 
in which the atoms are trapped in the radial dimensions 
by a coUimated laser beam, and trapped in the axial di- 
mension by another laser beam with a much larger waist. 
This combined optical potential would make independent 
adjustment of the trapping frequencies relatively simple. 
In addition, a pair of magnetic coils in Helmholtz con- 
figuration would provide the bias field needed to sweep 
the scattering length during the expansion without con- 
tributing to the trapping potential. 

By abruptly switching off the axial trapping poten- 
tial, the cloud would begin to expand freely along the 
axial dimension. The radial confinement would then be 
reduced by modulating the optical trap with a pair of 
crossed acousto-optic modulators, which can be used to 
produce a shallower, time-averaged potential which pro- 
vides the desired trapping frequency. With this tech- 
nique, we would expand the cloud until its optical depth 



is low enough, and its radial dimensions large enough, 
that atom noise correlations become measureable. At 
this point, the atoms would be imaged in situ using ab- 
sorption imaging. By appropriately choosing the optical 
transition and polarization for the imaging beam, atoms 
in the different spin states could also be distinguished 
and imaged independently in rapid succession. In this 
way, spatial correlations between atoms in the different 
spin states would be detected. 

We now make a few final remarks. The possible can- 
didate systems for this method would include the dilute 
single species Bose gas and the two spin or mass species 
Fermi gas at any value of the scattering length. In this 
latter case it is not necessary that the the numbers of 
atoms in each spin state be equal and so it could be also 
applied to the polarized Fermi gas. On the other hand it 
would be perhaps difficult to apply it to mixtures where 
more than one scattering length is present since we might 
lose the ability to control each of them simultaneously. 

We relie of course on the existence of an accessible 
Feshbach resonance and the application of a controlling 
magnetic field. Therefore, if the trapping potential is 
magnetic in origin then the two effects (trapping and de- 
tuning from resonance) must be handled together. Fi- 
nally we must keep in mind that the true potential seen 
by the gas includes the effect of gravity. When the gas 
begins to expand it will also begin to fall. This could 
introduce a distortion in the potential seen by the gas 
since it would move from its original position in the trap. 
To avoid this problem one possible solution would be to 
introduce a compensating magnetic field gradient. 

In conclusion, we began by discussing the importance 
and difficulty of measurements of correlation functions 
at short distances which can usually only be made af- 
ter expansion. We then show that, apart from the cases 
of the ideal and isotropic unitary gases, the second and 
higher order diagonal correlation functions do not ex- 
hibit scaling behaviour during the expansion of the gas 
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FIG. 1: Radial trapping frequency ujj_(t) from Eq. (|10|) for 
ujz(O) /lu±{0) = 0.5 (continuous line) and ijJz{0)/uj±{0) =0.1 
(dashed line) . Note that the t = discontinuity is appreciable 
in the first case (13%) but neghgible in the second (0.5%). 
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and so scaling cannot be used to relate the measured val- 
ues after expansion to those existing in the trap. We 
directly trace this violation of scaling to the existence of 
the Bethe-Peierls boundary conditions imposed on the 
many-body wave function due to the s-wave scattering 
between atoms. To address this problem we propose an 
experimental procedure that fully restores scaling to the 
diagonal correlation functions for systems originally in 
anisotropic traps and interacting via an arbitrary value 
of the scattering length which relies on a well-known 



ansatz for the many-body wave function. This ansatz 
is valid provided that the trapping potentials along dif- 
ferent directions can be changed independently and that 
the scattering length characterizing the interaction can 
be changed as a function of time with the use of a Fcsh- 
bach resonance. 
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